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Abstract 

Consider a relay cascade, i.e. a network where the source node, the sink node and a certain number of intermediate relay 
nodes are arranged in a hne. We assume that adjacent node pairs are connected by error-free (g + l)-ary pipes. The following 
communication scenario is treated. The source and a subset of the relays wish to communicate independent information to a 
. . . common sink under the condition that each relay in the cascade is half-duplex constrained. We introduce a simple channel model 

' for half-duplex constrained links and provide a coding scheme which transfers information by an information-dependent, non- 

' deterministic allocation of the transmission and reception slots of the relays. The coding scheme requires synchronization on the 

, symbol level through a shared clock. In the case of a relay cascade with a single source, the coding strategy is capacity achieving. 

^S) . Numerical values for the capacity of cascades of various lengths are provided, and it turns out that the capacities are significantly 

higher than the rates which are achievable with a deterministic time-sharing approach. If the cascade includes a source and a 
certain number of relays with their own information, the strategy achieves the cut-set bound when the rates of the relay sources 
fall below individual thresholds. Hence, a partial characterization of the boundary of the capacity region follows. For cascades 
composed of an infinite number of half-duplex constrained relays and a single source, we derive an explicit capacity expression. 
Remarkably, the capacity for g = 1 is equal to the logarithm of the golden ratio. We finally show that the proposed coding strategy 
is superior to network coding in the case of the wireless, half-duplex constrained butterfly network. 

Index Terms 



Half-duplex constraint, relay networks, network coding, timing, constrained coding, capacity, capacity region, method of types, 
golden ratio. 

I. Introduction 



C/3 

(N 

K*" ' A relay cascade is a network where the source node, the sink node and a certain number of intermediate relay nodes are 
arranged in a line. In this paper we consider the problem that a source node and an arbitrary but fixed number of relay 
ly-j nodes from the cascade wish to communicate independent messages to a common sink under the condition that each relay is 
^—1 half-duplex constrained, i.e. is not able to transmit and receive simultaneously. Throughout the paper, we assume that adjacent 
node pairs are connected by eiTor-free {q + l)-ary pipes. This approach allows us to gain a better understanding of half-duplex 
constrained transmission without having to distinguish which effects are due to channel noise and which result from the 
half-duplex constraint. Moreover, the problem becomes more feasible since combinatorial arguments can be used instead of 
statistical arguments. 

^ I How could we construct an efficient coding scheme which takes the half-duplex constraint into account? A first approach 
k>( ' would probably be to define a protocol such that the time-division schedule is determined a priori. Under this assumption, the 
rN capacity or rate region of various half-duplex constrained relay channels 111, O and networks |4| has been determined. We 
will show that time-sharing falls considerably short of the theoretical optimum or, conversely, higher rates are possible by an 
information-dependent, non-deterministic allocation of the transmission and reception slots of the relays. 

The meaning of information-dependent allocation scheme is illustrated in following example. Let Wq = {0, . . . , 7} be a 
message set. In each block i = 1,2, . . . of length 4, the source wishes to communicate a randomly chosen message wo{i) £ Wq 
to the destination via a single half-duplex constrained relay node. The alphabet of both source and relays equals {0, 1,N} 
where "N" indicates a channel use without transmission and {0, 1} is a g = 2-ary transmission alphabet. Let xo(i) be the 
codeword chosen by the source encoder to represent wo{i) in block i and let xi(i) indicate the codeword chosen by the relay 
encoder for representing wo{i — 1) in block i. The coding scheme is illustrated in Table I] The source encoder maps each 
message wo{i) to xo(i) by allocating the corresponding binary representation of wo{i), i.e. three bits, to four time slots. The 
precise allocation of the three bits to four slots is determined by the first two binary digits of codeword xo(i — 1). Based on 
the first two binary digits of the noiselessly received codeword xo(i), the relay encoder determines which time slot to use for 
transmission in Xi{i + 1). The binary value to be transmitted in xi(i + 1) is equal to the third bit in xo(i). Hence, the relay 
encodes a part of its information in the timing of the transmission symbols. Since the source encoder knows the scheme used 
to determine the relay's transmission slot, it can allocate its three new bits in xo(i + 1) to those slots in which the relay is able 
to listen. The sink then determines the message from the received relay codeword using both the position of the transmission 
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symbol and its value. In this example, a rate of 0.75 bit/use is achievable. By allowing arbitrarily long codewords, we will 
show that the strategy approaches 1.1389 bits per use which is also the capacity of the single relay cascade with half-duplex 
constraint when the transmission alphabet is binary. 



i? = 0.75 bit per use 


block i 


wo{i) 


xo(«) 


XI (j) 


jj)o(i) 


i = 1 


1 (001) 


OOIN 


NNNN 




i = 2 


2 (010) 


NOlO 


INNN 


1 


i = 3 


4 (100) 


INOO 


NONN 


2 


i = 4 


7 (111) 


llNl 


NNON 


4 













TABLE I 

The relay encodes a part of the information by the position of the transmission symbols. 



The example suggests that information encoding by means of timing is beneficial in the context of half-duplex constrained 
transmission. A similar example for q=l was shown in |5|. In |6|, Kramer applied the achievable decode and forward rates 
of the relay channel due to Cover and El Gamal [71 to a half-duplex constrained relay channel and noticed that higher rates 
are possible when the transmission and reception time slots of the relay are random. The randomness results from the fact 
that one can send information through the timing of operating modes. Timing is not a new idea in the information theoretical 
literature and has already been used in conjunction with queuing channels. Anantharam and Verdii showed JS] that encoding 
information into the distances of arrival to the queue achieves the capacity of the single server queue with exponential service 
distribution. The analog in discrete-time was analyzed in |9|. 

In Section we introduce a channel model which captures the half-duplex constraint in a simple way. A capacity achieving 
coding strategy based on allocating the transmission and reception time slots of a node relying on the node's previously received 
data is introduced in Section |III] The proposed strategy requires synchronization on the symbol level through a shared clock. In 
Section |IV] the performance of the coding strategy is analyzed yielding several capacity results. In the case of a relay cascade 
with a single source, it will be shown that the coding strategy is capacity achieving, i. e. approaches a rate equal to 

Cm-M = max min H{Y,\X,) (I.l) 

PXq.-.x^ l<i<m 

where m — 1 indicates the number of relays in the cascade and Xi and Yi are the sent and received symbol of relay i. If the 
cascade includes a source and a certain number of relays with their own information, the strategy achieves the cut-set bound 
given that the rates of the relay sources fall below individual thresholds. Hence, a partial characterization of the boundary of 
the capacity region follows. For cascades composed of an infinite number of half-duplex constrained relays, we show that the 
capacity is given by 

Coo(g) = log2 + bits per use. (1.2) 

Remarkably, Coo(l) is equal to the logarithm of the golden ratio. In Section [V] the capacity results are applied to various 
special cases. In particular, we transform (II.lt into a convex optimization program with linear objective and provide numerical 
solutions for Cm^i{q) for different values of m and q. In the case of a single relay channel with a source and a relay source, 
an explicit expression of the cut-set bound and of the achievable segment on the cut-set bound will be stated. We finally show 
that the proposed coding strategy is superior to network coding in the case of the wireless, half-duplex constrained butterfly 
network. 



II. Network Model 

We consider a discrete memoryless relay cascade as depicted in Fig. [T] The underlying topology corresponds to a directed 
path graph in which each node is labeled by a distinct number from V = {0, . . . , m} with m > 0. The integers and m 
belong to source and sink, respectively, while all remaining integers 1 to m — 1 represent half-duplex constrained relays, i. e. 
relays which cannot transmit and receive at the same time. The connectivity within the network is described by the set of 
edges £ — {{i, i + 1) : < z < to — 1}, i.e. the ordered pair {i, i + represents the communications link from node i to node 
i + 1. The output of the zth node, which is the input to channel + 1) is denoted as Xi and takes values on the alphabet 
X — {0, . . . ,q — 1}U {N} where Q — {0, . . . ,q — 1} denotes the g-ary transmission alphabet while "N" is meant to signify a 
channel use in which node i is not transmitting. The input of the ith node, which is the output of channel (i — 1, i) is denoted 
as Yi. Each message wq, sent via multiple hops from node to to at a transmission rate Rq, is uniformly drawn from the 
index set Wg — |l, 2, . . . , 2"^°} where n is the block length of the coding scheme. Besides forwarding previously received 
information, an arbitrary but fixed number of relay nodes also act as sources, i. e. each relay v E Vg intends to transmit its 
own messages at rate Ra(v) from 'Wa{v) — {I7 2, . . . , 2"^°<") } to the destination, where Vs summarizes all relays with their 
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Fig. 1 . A noiseless relay cascade (top) and the linic model illustrated by means of feedback. If relay i is transmitting, the switch is in position 1 otherwise 
in position 2. 

own messages and source node 0. Note that the bijection a : Vg ^ {0, . . . , |Vs| — 1} numbers the elements of Vs according 
to the order in which they appear in the cascade. Again, the transmission involves a multi-hop scheme since the information 
flow associated with a particular message Wq,(^) € yVa{v) has to pass all nodes between the corresponding (relay) source v 
and the destination. 

The output symbol of channel {i ~ is given by 



Y, 



if X,=N 
if X, e Q 



(11. 1) 



where 1 < i < m. Our channel model (III. Il l captures the half-duplex constraint as follows. Assume relay i is in transmission 
mode, i.e. Xi € Q. Then relay i hears itself (Yi = Xi) but cannot listen to relay i — 1 or, equivalently, relay i and relay i — 1 are 
disconnected. However, if relay i is not transmitting, i. e. Xi = N, it is able to listen to relay i — 1 via a noise-free {q + l)-ary 
pipe {Yi = Xi_i). The sink listens all the time, i. e. x^ is always equal to N, and therefore its input is given by Y„i = Xm-i- 
Another interpretation of the channel model is that the output of relay i controls the position of a switch which is placed at 
its input. If relay i is transmitting, the switch is in position 1 otherwise it is in position 2 (see Fig. [T]). Since a pair of nodes 
is either perfectly connected or disconnect, we obtain a deterministic network with p(yi, . . . , ym\xo: ■ ■ ■ , Xm) S {0, 1}. 



III. A Timing Code for Line Networks with Multiple Sources 
A. General Idea and Codebook Sizes 

A coding strategy is introduced which relies on the observation that information can be represented not only by the value of 
code symbols but also by timing the transmission and reception slots of the relay nodes. The strategy requires synchronization 
on the symbol level through a shared clock. 

The encoding technique applied at the source and the relays is as follows. The source uses a {q + l)-ary alphabet X 
for encoding without transmitting information in the timing of the symbols while relay i represents information by taking 
rii transmission symbols per block from the q-ary transmission alphabet Q combined with allocating the rii symbols to the 
transmission block. Then, at most " J different sequences x„i_i of length n can be generated by relay to — 1 where 

qUrr^-i (Jenotes the number of distinct sequences when the g-ary symbols are located at fixed slots while (^^ " ) denotes the 
number of possible slot allocations. Due to the half-duplex constraint, the effective codeword length of relay to — 2 reduces 
to n — n„i_i. This results from the fact that relay to — 1 cannot pay attention to relay to — 2 when relay m — 1 transmits and, 
therefore, the number of length n sequences producible by relay m — 2 is at most g"'"-^ ("^""^^). The same argument holds 
for each relay in the cascade. In general, relay i S {1, . . . , to — 1} is able to generate q"' (" "'^0 distinct sequences where 
rim = since the sink listens all the time. Finally, the effective length of the source codeword is n — ni which enables the 
source to generate {q + different sequences xq. 

Next, the maximum size of Wq, Wi, . . . , W|Vs|-i is given. Since the node with the smallest number of available sequences 
is obviously a bottleneck in the cascade from source to destination, we immediately obtain an upper bound on |Wo| which is 



iWol <min (q + 1) 



, min q ' 

l<i<m-l 



TLj 



(III.l) 



Both the source and the relay sources choose their messages uniformly and independent of each other Hence, relay source v 
is required to have nfio^ 1^*1 sequences available in order to represent an arbitrary set of arriving messages together with its 
own message. Consequently, for a]\ v ^ Vs \ {0} we obtain 



|Wa(.)| < 




mm q ' 

v<i<rn— 1 



n - n,+i 



(III.2) 



where the minimization in ( IIII.2I I yields the bottleneck, 
and the destination. 



e. the least number of available sequences, between relay source v 
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The transmission rate of (relay) source w e Vs is defined (in the standard way) as Ra(v) — log Assume 

the cascade contains Vg > 2 nodes with their own information. Further, assume that X^fi?"^ i?i > where Vs 9 f > 0, i.e. 
relay source v has to forward external information. Then the external information, which flows from node u — 1 to node v in 
block 6—1, has to determine the slot allocation used by node v in block b, i.e. the transmission pattern selected by node v is 
not allowed to depend on Wq(„)(6). Otherwise, node v — 1 would not know when transmission in block b is possible without 
collision where collision means that at least one transmission of two adjacent nodes occurs in the same time slot. Thus, the 
message sets have to satisfy for all v £ Vs \ {0} 

l>V.(.)l<| It^- . (III.3) 

Note that the bottom constraint is already contained in ( IIII.2I ). 

B. Example 

We now illustrate the ideas introduced in the previous section by constructing a code for a relay cascade with V = {0, . . . , 3} 
and Vs = {0, 2}. The transmission alphabet is binary, i. e. q = 2, and the code parameters are n = 4, ni = 1, n2 = 2. According 
to ( IIII.ll i to ( IIII.3I ), the maximum size of the message sets is {Wol = |Wi| = 4, which corresponds to a sum rate of 1 bit per 
use. Table HH depicts possible codebooks Co, Ci, C2 for nodes 0, 1 and 2, respectively. The codewords in the first row are used 
for representing source index 0, the codewords in the second row for representing source index 1 and so forth. The last row 
emphasizes that a codeword Xk {i) £ Ck, which is sent in block i by node k, represents a message wo(i — k) injected by the 
source encoder in block i — k with k £ {0,1,2}. 

Let us first consider €2- Four out of six possible transmission patterns are shown where a binary transmission slot is marked 
with B £ {0, 1}. Node 2 uses each transmission pattern for representing a particular source message wq £ Wq. Node 2's 
message wi £ Wi is encoded by the transmission symbols B. Each transmission pattern is identified with a unique color 
r £ {a, b, c, d}. 



wo 


Co 


Ci 


Ca 



1 
2 
3 


NONN e ONNN / ONNN g 
NINN e INNN / INNN g 
NNON e NNON / NONN g 
NNIN e NNIN / NINN g 


ONNN (a, e) NONN (6, /) ONNN (c, e) NONN {d, f) 
INNN (a, e) NINN (6, /) INNN (c, e) NINN (d, f ) 
NNON (a, g) NNNO {h, g) NNNO (c, g) NNON {d, g) 
NNIN (a, g) NNNl {h, g) NNNl (c, g) NNIN {d, g) 


NBNB a 
BNBN b 
NBBN c 
BNNB d 




wo{i) M xo(i) 

Example 


wo(i - 1) ^ xi(i) 
TABLE II 

CODEBOOKS FOR SOURCE, RELAY AND RELAY SOURCE. 


wo{i - 2) X2(«) 



Next, Ci is considered. Since node 1 knows the message wq to be forwarded by node 2 as well as codebook C2, it can always 
figure out both time slots ti and ^2 in which node 2 listens. Let xu^ , xu^ £ {0, 1, N} denote the symbols used by node 1 in ti 
and t2 for encoding a particular source message. The following mapping wq i— > {xu-^ , 2:1*2) is chosen: 1-^ (0, N), 1 1-^ (li N), 
2 (N, 0), 3 1-^ (N, 1). By allocating each of the four values of {xit^ , xu^) to the listen slots of pattern r £ {a, b, c, d] and, 
further, by requiring that node 1 is quiet when node 2 sends a binary symbol B, we obtain the codewords in Ci which are 
colored by (r, s). Color s £ {e,f,g} labels the resulting transmission patterns in Ci. Color r helps node 1 to pick the new 
codeword from the correct column, i.e. when node 2 uses a pattern with color r node 1 uses in the same transmission block 
a codeword whose first color is equivalent to r. 

Finally, we consider Cq. In each transmission block, source node can use three time slots ti, t2 and for encoding since 
node 1 sends once per block. Let xotj , xota , xotj £ {0, 1,N} denote the symbols used by node for encoding a particular 
message wq £ Wq. We use a similar mapping as before, i. e. wq t-^ (a^oti 1 2^0*2 1 ^ots ) is given by 1-^ (0, N, N), 1 1-^ (1, N, N), 
2 1-^ (N,0,N), 3 I— > (N, 1,N). Now, by allocating all possible values of (xoti , 20*2 7 2^0*3) to the listen slots of codewords in 
Ci whose second color is s and, further, by requiring that node is quiet when node 1 transmits, we obtain all codewords 
in Co which are colored with s. It should be noted that merely four from 27 possible sequences are used in the mapping 
Wq (xoti , 2:0*2 , 2:0*3). Hence, Co could be designed such that node is able to send [27/4J additional messages to a sink 
at node 1 at a rate of 0.6462 bit per use. In summary, the source encoder applies the following strategy. Based on message 
'Wo{i — 2), the first color r of codeword xi(i) is determined. Subsequently, based on this information, the source determines 
the second color s of Xi(i) by means of 'Wo{i — 1). This color tells node from which column the new codeword has to be 
picked, namely from a column whose codewords are colored with s. The precise choice within the picked column depends on 
the new source message wo(*)- 
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C. Rate Region 

We now turn towards gaining an achievable rate region TZ from the expressions derived in section IIII-AI The following 
abbreviations are used for the portion of time in which relay i listens or transmits: pi = n~^{n — rii),0 < rn < n, and 
pi = 1 — Pi where pi — px^ (N), 1 < i < m. 

The method of types 1 10 1 provides important tools for relating combinatorial expressions to information theoretic expressions. 
An example very useful for the problem considered here is 1,11. Th. 1.4.5] 



n log 



H {pi) + o(l) for n — !■ oo 



(111.4) 



where H{pi) denotes the binary entropy function evaluated at pi — n ^{n ~ Ui). By (IIII.4l i and Ri = log |yVi|/n, we obtain 
from (llll.lb to (IIII.3b for n ^ oo 



Ro < mill \ pi \og{q + 1), min (pi log q + Pi+iH (piP^+i) + o(l) ) } 

I l<2<m— 1 \ J 



a{v) 



Vi?, < min {p,logg+p.,+ii/(KP,+i) +o(l)}, Vv G V, \ {0} 

' ^ 'i?<2<m— 1 



2^0 



a{v)-~-l 



Rc.{v) < Pv\ogq, if ^»>0' Vz;GV,\{0}. 



(III.5) 
(III.6) 

(111.7) 



Further simplifications are possible by taking into account the optimal structure of the marginal distributions pxoXi , ■ ■ ■ , 
PXm~iX„^ as shown in Tables |III] and |IV] The zero probabilities in Table |III] and |IV] result from following consideration. 



^0 



9-1 



N 



— —■ + 

■■■ : 

pi/q ■■■ Pi/q pi/Cq + i) 

TABLE III 
OPTIMAL PXqXi- 



i-1 



X, 



9-1 



N 

Pi-l/9 



9-1 ■■■ Vi-i/l 

N pi/q ■■■ pi/q pi-pi-i 

TABLE IV 

Optimal pXi_iXi for 2 <i<m. Note THATpm = 1. 



Assume relay i is transmitting, i. e. Xi E Q. According to the underlying channel model, relay i is not able to detect the input 
of node i — 1 and, consequently, node i — 1 should not transmit when node i transmits. Or, to be more precise, a channel input 
pair {xi^i, Xi) is negligible if it produces the same channel output pair (j/i, j/i+i) as another channel input pair and this with 
the same probabilities. Hence only one non-zero entry remains in each of the first q columns of Table Hill and IIVI whereas the 
assignment of the non-zero entry within a column is not unique from an information theoretic viewpoint. However, from an 
engineering point of view, the assignment as depicted in both tables is reasonable since node i — 1 should not transmit and, 
therefore, waste transmit power when its input cannot be detected by node i. As a simple consequence of the zero probability 
assignment, we have the relation pXi^iXii^, k) — PXiXi+i{k,N) for all fc G Q and 1 < i < m — 1. 

Let us now address the remaining values in Table Hill and IIVI First, consider the time slots in which the first relay listens. 
During this fraction of time, the source should make optimum use of the channel by encoding with uniformly distributed input 
symbols. Hence, pxoXi{k,N) = pxgXiil,^) for all fc,Z G X. By taking the relative frequency of the transmission symbols 
into account, we have pxoXi{k,N) = Pi/{q + 1). Moreover, in order to achieve the maximum information flow from relay 
i — 1 to node i or, likewise, from the fact that a permutation of the transmission symbols Xi-i G Q obviously yields the same 
information flow, we can choose pXi-iXi {k, N) = PXi-iXt {I, N) for all k,l G Q. Due to the relative frequency of transmission 
symbols within a block, we have pXi-iXi {k, N) = Pi-i/q for all /c G Q where 2 < i < m. 
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It is now fairly easy to check that the following equalities hold 

HiYi\Xi) = pilog(g + l) (III.8) 

H{Y,+i\X,+i) = p,logq+p,+iH {p,pr^,) , l<z<m-l. (III.9) 

From (IIII.5I ) to jlll.TI l together with dlll.SI ) and ( IIII.9b . we obtain 

7^ Co IJ 7^p^ (III.IO) 

where 7?.p for n ^ oo is given b}Q 



I p:=PA'o...x„ 



7^p= <jR:yi?fc< min V« G V, 

' k=0 

( q(u)-1 

f|iR:i?,(„) <p,„loggV ^ i?, = 0, Vi- G V, \ {0} 



(III. 11) 



and Co(— ) denotes the convex hull, i.e. takes timesharing into account. Note that R is a Vs -dimensional rate vector with Ri 
as its ith entry. 

IV. Capacity Results 
In this section we shall investigate the optimality of the coding strategy. 

Theorem I: A part C of the capacity region C of a noise-free relay cascade with |Vs| sources and m — 1 half-duplex 
constrained relays is given by 

C = M {ll:yRk< min H{Y,\X,), G V, i (IV.l) 

^ I v + l<i<m I 



k=0 



i.e. if the elements of R satisfy 



a{v)-l 

R : Rciv) < p. log g V y i?, = 0, Vi; G V, \ {0} I (IV.2) 

i=0 J 

then (II V. Il l yields the corresponding boundary points of C for some joint distribution pxo...Xm- 

Proof: We first show which subset of the network cuts is sufficient for the considered line network. Recall that an upper 



bound on the sum rate J2k=o is given by LI 2. chap. 14.10] 



Vi?fe< max mmIiX,,Xs^:Ys,Y,n\Xs), (IV.3) 

^ — ' Px„...x,„ seM 

fc=0 

where M. = V {{v + 1, . . . , m — 1}) and S'^ is the complement of in {u + 1, . . . , m — 1}. Since our network model is 
deterministic, (IIV.3I ) simplifies to 

a(v) 

y^Rk< max imn H{Ys,Y^\Xs). (IV.4) 

k=0 

Now assume that S is nonempty and let i G {w + 1, . . . , m — 1} denote the smallest integer in S. By the chain rule for entropy, 
we can expand H{Ys,Ym\Xs) as 

HiYs,Y„,\Xs) = H{Y,\Xs) + H{Ys\[,}\Xs,Y,) + H{Y„,\Xs,Ys) (IV.5) 
> H{Y,\Xs). 

For each cut S with smallest entry i, a cut here called Si can be found such that H{Yg^,Ym\Xs.) is less than or equal to 
H(Ys, Ym\Xg). Simply choose Si :— {i, . . . ,m — 1}. This eliminates the second and third term on the rhs of ( |IV.5t due to 
the underlying channel model. Further, since S C Si we have H{Yi\Xg) > H{Yi\Xsi) and, thus, each non-empty cut S with 

'Note that the Landau symbols are neglected for the sake of simple notation. 
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smallest element i is dominated by Si in terms of delivering a smaller entropy value. Finally, S* = has to be considered in 
(|IV4] |: S* = yieldfl H{Y„,). To sum up, J2k=a is upper bounded by 

Rk < max min H{Y,\Xs,) (IV.6) 

fc=0 V m 

< max min H{Y,\X,) (IV.7) 

Px^...Xjn v+l<i<m 

where the last inequality follows from the fact that conditioning reduces entropy. Therefore, the cut-set bound C is given by 

( c,{v) ^ 

C= 11 iR -.y" Rk< min H{Y\X,), Vw G > . (IV.8) 

^ I ti+l<i<m I 

PXo...x„ k=0 ) 

Let C denote the capacity region, i.e. C C C. If we focus on rate vectors R whose elements satisfy iTV.li . then TZ' defined as 

7^'= U Up (IV.9) 

p:=PXo...x„ 

(i.e. TZ from (IIII.lOl l and (llll.llb without the timesharing points) equals C. Thus, TZ' — C under constraint ( IIV.2l i. ■ 

Corollary I: The capacity of a noise-free relay cascade with a single source-destination pair and m— 1 half-duplex constrained 
relays is given by 

C„_i(g) = max min H{Y,\X,) (IV. 10) 

PXo...x,„ l<i<m 

where q equals the number of transmission symbols. 

The capacity of a single source line network with an infinite number of half-duplex constrained relays is stated in Theorem|2] 

Theorem 2: For to — > oo, i. e. for an unbounded number of relays, and q transmission symbols, the capacity of the noise-free 
and half-duplex constrained relay cascade with a single source-destination pair is given by 



Cooiq) = log2 ( ^ + ^^'i + ^ ) bits per use. (IV. 11) 



Proof: Theorem |2] is proved in the Appendix. ■ 
Remarks: 

i) In order to achieve Coo (q) it follows from equation iA.6\ that each relay has to transmit t percentage of the time where 

T = 50(l F^= I . (IV.12) 

ii) ^00(1) = 0.6942 bit per use is equal to the logarithm of the golden ratio. Also remarkable. Coo (2) is exactly 1 bit per 
use. 



iii) The maximum achievable rates with time-sharing and, thus, no timing, are given by Rts{q) — log2 \/q + 1 bits per use. 
For q = 1,2 we have 0.5 and 0.7925 b/u, respectively. Since Coo{q) is obviously a lower bound on the capacity of 
each finite length cascade, a comparison of the time-sharing rates with Coo(l) and Coo (2) shows that time-sharing falls 
considerably short of the theoretical achievability for small transmission alphabets. For very large transmission alphabets 
the gap between the rates due to time-sharing and timing becomes negligible, i.e. lim^^oo (Coo (9) — Rts{q)) — 0. 

V. Numerical Examples 

In this section we shall provide numerical capacity results for various scenarios by means of Theorem[T]and Corollary [1] In 
particular, we show how to obtain the capacity of a half-duplex constrained relay cascade with one source-destination pair for 
an arbitrary number of relays. Further, in case of a single relay cascade with source and relay source, an explicit expression 
of the region due to Theorem [T] is derived. Throughout the section, the base of the logarithm is assumed to be two. 

2Note that H{Ym) = H(Yrr,\X„,). 
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Fig. 2. Graphical solution of optimization problem )V. It . 
A. One Source 

Let us first consider a relay cascade with V = {0, 1, 2}, = {0} and q ~ 2, i. e. source node intends to communicate 
with sink node 2 via the half-duplex constrained relay 1. By Corollary [T] and the optimum input pmf stated in Table HIH we 
have 

Ci(2)= max min {p^, (N) log2 3, . (V.l) 

PX0X1X2 

Problem jV.ll ) exhibits a single degree of freedom. Since the maximum does not occur in the maximum of one of the two 
(concave) functions, the problem is readily solved by finding a px^i^) which satisfies pxi(N) logj 3 = H{Xi) (see Fig. |2l). 
The optimum value for pxi(N) equals 0.7185 which yields 

Ci(2) = 1.1389 bits per use. (V.2) 

Remarks: 

i) Assume the relay does not have the capability to decide whether the source has transmitted or not, i. e. pxoXi (N, N) = 0. 
In this case an identical approach shows that the capacity equals 0.8295 bit per use, which is still greater than the time- 
sharing rate of log2 \/3 bit per use. 

ii) For g = 1, the outHned procedure yields Ci(l) = 0.7729 bit per use achieved by px^ (N) = 0.7729. The capacity value of 
this specific case has also been obtained in ITSl . Therein, the focus was not on half-duplex constrained transmission but 
on finding the capacity of certain classes of deterministic relay channels. In |5 |, the same channel model was considered 
and the author noticed that the capacity is greater than 0.5 bit per use though a half-duplex constrained relay is modeled. 
A simple coding scheme was outlined which approaches 2/3 bit per use. 

In order to compute Cm-i{q) for to > 2, we transform (llV.lOb into a convex program with linear cost function H{Yi\Xi) 
and convex equality constraints H{Yi\Xi) — H{Yi+i\Xi+i) = for all i £ {1, . . . ,to — IjH. The resulting program reads as 

maximize pi \og{q + 1) 

subject to pi \og{q + l)+pi log — h (pi+i - Pi) log — — = (V.3) 

qPi+i Pi+i 

l-^p,<0 
Pr e [0, 1] 

By adopting a standard algorithm for constrained optimization problems, the capacity Cm-i{q) was computed for various 
values of m. A brief summary is given in Table [V] 

^See proof in the Appendix why the constraints are satisfied with equality. 
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0.6954 b/u 


1.0020 b/u 


41 


0.6946 b/u 


1.0006 b/u 


101 


0.6943 b/u 


1.0001 b/u 


oo 


0.6942 b/u 


1 b/u 


TS 


0.5 b/u 


0.7925 b/u 



TABLE V 

Capacity results for various cascades composed of m - 1 half-duplex constrained relays. The transmission alphabet is either 

UNARY OR BINARY. ROW "TS" SHOWS THE CORRESPONDING TIME-SHARING RATES . 



B. Two Sources 

The considered relay network is characterized by V = {0, 1, 2}, Vs = {0, 1} and g = 2. By Theorem [U a part C of the 
capacity region is given b}0 

i?o < H{Xo\X^) (V.4) 
i?o + i?i < H{Xi) (V.5) 
Ri < Pi, ifi?o>0 (V.6) 



for some pxoXi- We will first derive an explicit expression for the boundary of C which is characterized by ( IV.4b and dV.SI) . 
Subsequently, C and TZ are given. 

Two cases have to be considered depending on whether an optimum input pmf for the source or the relay source is used. 
An optimum input pmf for the relay source is shown in Table [VTl It yields a maximum sum rate H{Xi) of logj 3 bits per use 
for all valid y (i.e. y e [0, 1/6]). When y varies from to 1/6 all points on i?i = log2 3 — i?o for < i?o ^ ■|log23 result. 



^0 







1 



N 




1 
N 






1/3 1/3 

TABLE VI 

Optimal relay pxqX-i which yields a sum rate of logj 3 bits per use. 



1/3 - 2y 



It remains to focus on the interval |log23 < i?o < 1.1389 bits per use. Under consideration of the optimum input pmf for 
source node (Table Ullb and ( IIII.SI l. we can express i?i — H{Xi) — Rq as shown in the second line of ( IV.7l i. Hence, the 
boundary of C is given by 

j log23-i?o, 0<i?o<|log23 
""'"l ^(Tif3) + (l"IS^)-^o, ilog2 3 <i?o< 1.1389 ■ (^-^^ 

In order to determine C, (IV.6I 1 has to be taken into account. For Rq > (IV.6I 1 yields an upper bound on Ri or, equivalently, 
a lower bound on Rq. This lower bound is given by the right hand side of 

H{Xo\Xi)>H{Xi)-Pi. (V.8) 

Since H{Xo\Xi) is linear in pi while H{Xi)—pi is concave in pi, the smallest value for the lower bound follows by finding a 
pi which achieves equality in ( IV. St . We obtain pi = 0.6091 what gives Rq > 0.9654 and Ri < 0.3909 bit per use. If Rq = 0, 
(IV.6I 1 is not valid anymore and we have i?i < log2 3. Thus, the boundary of C is given by ( IV.7l i for 0.9654 < Rq < 1.1389 bits 
per use together with the rate vector R = (0, log2 3) bits per use. TZ follows from C by taking the convex hull. In particular, all 
points on the connecting line between (0, log2 3) and (0.9654, 0.3909) bits per use are added. The three regions are depicted 
in Fig. [51 

The derivation reveals the following interesting fact. Even when the source transmits at a rate beyond the time-sharing rate 
of log2 Vs bit per use, the relay is still able to send its own information at a non-zero rate. 



'^H(Xi) is not considered on the right hand side of )V.4t since it already appears in jV.St 
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Fig. 3. C, C and Ti are depicted. Tlie points marked witli a star and a circle coirespond to (0, log2 3) and (0.9654, 0.3909) bits per use. C equals the 
region bounded by the solid curve. C is the region bounded by the dotted line and the fraction of the solid curve to the right of (0.9654, 0.3909) together 
with the isolated point (0, log2 3). H is the union of C' and the region bounded by the dash dotted line. 

VI. Extension to other Networks 

Relay cascades are fundamental building blocks in communication networks. Therefore, the results derived in the previous 
sections may be instrumental in order to determine the capacity of half-duplex constrained networks with more elaborate 
topologies. 

A. Wireless Trees 

Consider, for instance, the tree structured network depicted in Fig. |4] The root (node 1) wants to multicast information to 
all leaves (nodes 2 to 8) via four half-duplex constrained relays. We assume noise-free bit pipes (i. e. q ^ 1) and broadcast 
behavior at nodes with more than one outgoing arrow. The multicast capacity is limited by the capacity of the longest path 
in the tree which goes from node 1 to node 7 or 8. Hence, the multicast capacity in the considered example is equal to the 
capacity of a cascade containing two intermediate relay nodes, i. e. C2(l) = 0.7324 bit per use (see Table IVTl. 



1 




Fig. 4. A wireless binary tree. The multicast capacity is equal to C2(l) = 0.7324 bit per use. 



B. Wireless Butterfly 

Another example for a wireless butterfly network lfT4l is shown in Fig. |5] Nodes 1 and 2 intend to multicast information to 
sink nodes 4 and 5 via both a direct link and a half-duplex constrained relay node 3. Like before, broadcast transmission and 
bit pipes are assumed. All nodes with two incoming arrows behave according to a collision model, i. e. received information is 
erased if there was a transmission on both incoming links. By means of network coding (NC), 2/3 bit per use are achievable at 
the sink nodes, as is illustrated in Fig.|5](a). The (well-known) strategy is to send in the first time slot a binary symbol ui via 
broadcast a to nodes 3 and 4, in the second time slot a binary symbol U2 via broadcast b to nodes 3 and 5 and, subsequently, 
in the third time slot ui U2 via broadcast c to both sinks. However, under the usage of timing, at least 0.7729 bit per 
use is achievable as is illustrated in Fig. |5] (b). This results from the fact that information originating from node 1 can be 
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sent by means of timing at a rate of Ci(l) = 0.7729 bit per use (see Table |V]l concurrently on paths 1, (1,3), 3, (3, 4), 4 
and 1, (1,3), 3, (3, 5), 5. Equivalently, information originating from node 2 can also be sent by means of timing at a rate of 
Ci(l) = 0.7729 bit per use concurrently on paths 2, (2, 3), 3, (3, 4), 4 and 2, (2, 3), 3, (3, 5), 5. Hence, time-sharing of both 
source nodes yields a multicast rate of 0.7729 bit per use. It should be noted that the direct Unks (1,4) and (2,5) are not 
necessary in order to achieve Ri — R2 — 0.7729 bit per use, which suggests that the multicast capacity is even larger. 




Fig. 5. The wireless binary butterfly network. With network coding 2/3 bit per use are achievable. Timing yields Ci(l) = 0.7729 bit per use. 



VII. Conclusion 

The half-duplex constraint is a property common to many wireless networks. In order to overcome the half-duplex constraint, 
practical transmission protocols deterministically split the time of each network node into transmission and reception periods. 
However, this is not optimum from an information theoretic point of view, as is demonstrated in this paper by means of 
noise-free relay cascades of various lengths with one or multiple sources. We show that significant rate gains are possible 
when information is represented by an information-dependent, non-deterministic allocation of the transmission and reception 
slots of the relays. Moreover, we provide a coding strategy which realizes this idea and, based on the asymptotic behavior of 
the strategy, we establish capacity expressions for three different scenarios. These results may be instrumental in deriving the 
capacity of half-duplex constrained networks with a more elaborate topology. 



Appendix 

Proof of Theorem^ It is first shown that WpXi^iXi PXiXi+i\\oo ^ if to ^ 00 for all 2 < z < m— 1. The capacity 
series (Cm('z))meN is bounded (e. g. by and Ci{q)) and monotonically decreasing what follows from the fact that each new 
relay causes an additional constraint in the corresponding convex program ( IV.3l l. Hence, {Cm{(l))me'M is convergent, i. e. for 
every e > there exists an G N such that 

\Cra{q) - Cra+l[q)\ < e (A.l) 

for all 771 > A^. 

Further, the capacity achieving input pmf pxo...x^ in (II V. 101 ) yields H{Yi\Xi) = H{Yj\Xj) for all I < i, j < m. Assume 
this is wrong, i. e. there exist indices i,j such that H{Yi\Xi) > H{Yj\Xj). However, H{Yi\Xi), i.e. the transmission rate 
from node i — 1 to node i, can be decreased without forcing any of the remaining nodes to decrease their transmission rates. 
Simply reduce the fraction of a block used by node i — 1 for encoding until the desired rate is achieved and fill the remaining 
slots of the block with quiet symbols N. 

Hence, assuming the capacity achieving input pmf, we have for instance Cm{q) = H{X,n) and Cm+i{q) = H{X„i^i). 
Then, by dAlT ) 

|i7(X,„)-i/(X„+i)| <e (A.2) 

for all TO > A^ where 

1 — Pk 

H{Xk) ^ -Pklogpk - {I ~ Pk) ^og — - — , fc = TO,TO+l (A.3) 

(see Fig. |2]for a plot of H{Xk) for q — 2). Two cases can appear in ( IA.2b when e approaches zero: pm and Pm+i are arbitrarily 
close to a common point p or pm and p„i+i are arbitrarily close to two distinct points p' , p" . 

We note that if the second case occurs, p' + p" is not allowed to be smaller than one since otherwise negative probability 
masses would result (see Table II VI ). However, p' +p" is always smaller than one what can be seen as follows. First, note that the 
maximum of H{Xk) is at l/{q+l). Hence, without restriction we can assume that p' < 0.5 and p" > 0.5 (otherwise p'+p" < 1 
a priori). Since the first derivative of H{Xk) is point symmetric with respect to (0.5, — logg), we have 0.5 — p' > p" — 0.5 
what yields p'+p" < 1. 
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Hence, only the first case is valid, i. e. \pra — Pm+i\ ^ as m — > cx). But this implies WpXi^iXt — PXiXi+i\\oo for all 
i e {2, . . . ,m - 1} and, thus, \H{X,^i\X,) - H{Xj^i\Xj)\ -> for all 2 < i, j < m - 1 as m -> oo. 
In the final step of the proof, we show that the capacity Coo (?) is equal to the maximum in p of 

= -(1 - p) log - {2p ~ 1) log (A.4) 

qp p 

where i > 1. Elementary calculus yields 



max H{X,^,\X,) = log I ^ + '^^'^ + ^ ) (A.5) 



achieved at 



It remains to show that H{Xo\Xi) evaluated at p*, i. e. 

H{Xo\X,)\p. = ^ (l + ^1=^) log('? + 1) (A.7) 
is always greater or equal to (IA.5I) . This is satisfied if 

(, + lfO"^)>^ + ^ (A.8) 

or, more strictly, 

(, + l)K^^^)>i±M±I. (A.9) 

Using the substitution 

2V9TT 



g= ^ J-^ (A- 10) 



in ( IA.9I 1 yields 

(2g)-''>g+l. (A. 11) 

jA.l 11 1 is satisfied for all q S [0,0.2] what can be seen as follows. First note that ( lA.l lb is satisfied for q = and q = 0.2. 
Since {2q)~'' is concave due to a non-positive second derivative in the considered domain, dA.llb is valid for all q e [0, 0.2]. 
Thus, ( IA.8I ) is true for all q > 5. The vahdity of ( IA.8I ) for the remaining q £ {1, ... ,5} can be easily checked. ■ 
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